Quantum dots in magnetic fields: Phase diagram and broken symmetry of the 

Chamon-Wen edge 
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Quantum dots in magnetic fields are studied within the current spin density functional formal- 
ism avoiding any spatial symmetry restrictions of the solutions. We find that the maximum density 
droplet reconstructs into states with broken internal symmetry: The Chamon-Wen edge co-exists 
with a modulation of the charge density along the edge. The phase boundaries between the polar- 
ization transition, the maximum density droplet and its reconstruction are in agreement with recent 
experimental results. 
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Quantum dots are small electron islands made by lat- 
erally confining the two-dimensional electron gas in a 
semiconductor heterostructure. Such nano-sized systems 
attracted much interest since the techniques in their fab- 
rication developed far beyond mesoscopic dimensions [Q . 
Vertical quantum dots can nowadays be made so small 
that they show atom- like behavior |g] : shell structure and 
Hund's rules determine the electronic properties. 

Much experimental effort focussed on systematically 
mapping the magnetic field dependence of the chemical 
potential obtained from single-electron capacitance spec- 
troscopy ^. As a bias is appHed to the gates, single 
electrons tunnel into the quantum dot when its chemi- 
cal potential IJ.{N, B) (which depends on the number of 
confined electrons N and the magnetic field strength B) 
equals the Fermi energy in one gate electrode. First ex- 
periments along these lines were performed by Ashoori et 
al. [|j and later Klein et al. Q. Recently, Oosterkamp 
et al. js) systematically extended the measurements to 
stronger fields B and larger sizes N . Cusps and steps in 
yLt(iV, B) were found to clearly separate different ranges 
of magnetic fields. From a comparison to results of ex- 
act diagonalization studies 1^ these patterns were identi- 
fied with phase transitions in the droplet: they occur at 
magnetic fields for which the ground-state charge distri- 
bution of the dot changes, defining sharp phase bound- 
aries. The points at which a complete polarization of 
the electrons occurs mark the beginning of the so-called 
Maximum Density Droplet (MDD) phase. This new state 
suggested by McDonald, Yang and Johnson Q is a ho- 
mogeneous droplet in which the density is approximately 
constant at the maximum value po — (Stt^^)^^ that can 
be reached in the lowest Landau level. (£b = yJh/eB 
is the magnetic length). In the spin-polarized MDD the 
electrons occupy adjacent orbitals with consecutive angu- 
lar momentum. This compact occupation of states max- 
imizes the electron density. The stability of the MDD 
is determined by a competition between the kinetic and 
external confinement contributions to the total energy, 
and the Coulomb repulsion of the electrons. The former 
would favor the MDD structure up to infinite fields: with 
increasing _B, the droplet would decrease in radius such 
that close to the dot center it could maintain a density 
corresponding to filling factor one in the bulk limit [Q. 



This, however, is inhibited by Coulomb repulsion: At a 
sharply defined transition point, the charge density dis- 
tribution of the droplet reconstructs Chamon and 
Wen Q found from Hartree-Fock calculations that at the 
edge of a quantum Hall liquid with bulk filling factor one, 
a stripe or ring of electron density breaks off from the ho- 
mogeneous bulk or the MDD phase. 

In this Letter, we show that the so-called Chamon-Wen 
edge actually is a ring of nearly localized electrons sur- 
rounding the MDD (see Fig. || below). Moreover, we find 
that the calculated phase boundaries are in good quali- 
tative agreement with recent experimental results ||5[| . 

As exact diagonalization techniques |^ are limited to 
small particle numbers, mean-field methods are needed 
to simplify the complicated many-body problem. We ap- 
ply the so-called current spin density functional theory 
(CSDFT) including gauge fields in the energy func- 
tional. In contrast to the HF calculations by Chamon 
and Wen §] or recent applications of CSDFT we 
avoid any spatial symmetry restrictions of the mean field 
solution. 

As a basic model for a quantum dot one usually con- 
siders N interacting electrons of effective mass m* con- 
fined in a two-dimensional harmonic trap. A homoge- 
neous magnetic field B = Bgz is applied perpendicular 
to the x-y-plane in which the electrons are confined by 
the external potential V = m*uP'r'^ jl. In a symmetric 
gauge, the external vector potential acting on the elec- 
trons is A = B/2{-y,x). For the details of CSDFT we 
refer to the work of Vignale and Rasolt and recall 
here only the most important steps. The self-consistent 
Kohn-Sham type equations in CSDFT read 



2ra* 



2m' 



(p • ^ + A-Tp) + V5 



(We have dropped the arguments r for simplicity). The 
index i labels the eigenstates with spin 5 — (t, [). We de- 
fined effective vector and scalar potentials A:= A-\- Axe 
and Vs := {e'^ /2m*)A'^ + Vs + Vh + V^cS- Here, Vh is the 
ordinary Hartree potential and Vs = V + {—)g*^BB/2 
is the external potential, including the Zeeman energy, 
(/is = eh/{2me) is the Bohr magneton). The exchange- 
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correlation vector and scalar potentials are 

^ 1 / <9 a[pexc(p^,7)] d d[pey,c{p5n)] \ 

p \ dy 97 ' dx dj / 

and 

_ d[pe^cips,l)] e 
9P5 P 

where p is the particle density p = p-\ + Pi with = 
J2i l^i<5p • The paramagnetic current density is given by 
jp = -i?i/(2TO*)Ei5[«'*5V4'.5 - *»5V**5] and the real 
current density equals j = jp + {e/m*)Ap. In essence 
and as applied in CSDFT, the exchange-correlation en- 
ergy depends on the so-called vorticity 7 = V x (jp/ p)\^ 
of the wave function. In the bulk the total current den- 
sity must be zero, and thus 7 = —eB/m*. It is this 
relation that allows to use the interpolation formulae for 
the exchange energy per particle Cxc for the homogeneous 
bulk in CSDFT by replacing B {m*/e)\-f\. Mak- 
ing use of the local spin density approximation (LSDA), 
the exchange energy per particle Cxc is parametrized in 
terms of the local particle density p, the spin polarization 
^ = (p| — Pi)/p and the filling factor 1/ = 2TThp/ni*\j\. 
We used the expression exc(p, C;^^) — ^^c{p)^~^^^'' + 
ejf (p,0(l - e-^('')), where f {v) = l.bv + 7i/^ This 
form interpolates between the infinite magnetic field 
limit e^(p) = — 0.782-y27rpe^/47reoe and zero field limit 
e^^{p, \), for which we use the Tanatar-Ceperly func- 
tional and generalize it to intermediate polarizations [ p5[ . 
For V < 0.9, the interpolation in e^cip, ^) follows 
closely the results of Fano and Ortolani |lq| for polar- 
ized electrons in the lowest Landau level, and saturates 
quickly to the zero field result for v > 1. 

The LSDA has been shown to describe very accurately 
addition spectra [^-^ in weak or zero fields (110 ■ Here 
we concentrate on the polarization transition and be- 
yond. We use a plane wave basis to solve the single- 
particle Kohn-Sham equations self-consistently. The 
practical computational techniques which we found nec- 
essary to obtain converged solutions of the CSDFT mean 
field equations are described in Rcf. ||l8|. For the mate- 
rial parameters we choose the typical GaAs values for 
the effective mass m* = 0.067me, the dielectric constant 
e = 12.4 and the reduced Lande g-factor g* — 0.44. 
(This yields an effective Bohr radius of a*g = 9.79nm.) 
The strength of the external confinement is set to Tioj = 
4.192 iV~^/'*meV . For different TV, this convention keeps 
the average electron density in the droplet approximately 
constant, in this case roughly corresponding to a (two- 
dimensional) Wigner-Seitz radius = 2a|j. 

As a first example we study iV = 20 electrons confined 
in the above specified GaAs dot. At a field of about 2.3T, 
only one minority spin is left in the dot center. The 
droplet becomes completely polarized for fields larger 
than about 2.4T. Shortly after the polarization point the 
MDD is formed. Its density profile has perfect azimuthal 
symmetry and is nearly constant inside the droplet corre- 
sponding to filling factor one in the bulk: all electrons are 



in the lowest Landau level and all single-particle states 
are occupied by one spin down electron with successive 
angular momenta to = 0, l,...,iV — 1 M]. The total or- 
bital angular momentum then is M = ^N{N — 1), which 
yields M =160. As in the inner regions of the MDD 
the filling factor equals one, the density rises with in- 
creasing magnetic field, and the whole droplet decreases 
in radius. For N =20 the MDD is stable for fields up 
to B ^2.9T, where the reconstruction of the charge den- 
sity distribution begins. The density distribution of the 
MDD at i? = 2.9T is shown in the upper panel of Fig. |^. 



N=20 




FIG. 1. Left: Self-consistent densities for a 20-electron 
GaAs quantum dot shortly before edge reconstruction 
at B = 2.9T (top), forming the broken-symmetry Cha- 
mon-Wen edge at 3.0T (middle) and being fully recon- 
structed at 3.5T (bottom). Right: Angular momentum 
occupancies Pm as defined in the text 

We can gain information about the angular momentum 
occupancy in the droplet by projecting the Kohn-Sham 
single-particle wave functions l^'ij) on Fock-Darwin l2C| ] 
states \nrn) with radial quantum number n and good an- 
gular momentum to. Summing over all occupied states 
i < N, we obtain the total angular momentum occu- 
pancy P„i = Z^i.n ■^itli * ^ ^- (In strong 
fields, only n — states give a significant contribution.) 
The values of Pm are plotted with bars for each m- value 
in the right column of Fig. |l|. As expected, for the MDD 
all angular momenta up to — 1 have occupancy one. 

Increasing the magnetic field effectively compresses the 
electron states, such that states with constant m decrease 
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in radius. Consequently, the density of the MDD would 
continuously increase with B. This, however, is inhib- 
ited by the Coulomb interaction. The droplet thus re- 
distributes 1^ its density over the dot area, taking ad- 
vantage of moving electrons from lower to higher angular 
momentum states and leaving some partly unoccupied 
states for m < N ~ 1. Thus, M > iV(V - l)/2 and the 
Chamon-Wen edge is formed . The corresponding den- 
sity and angular momenta at B =3.0T are shown in the 
middle panel of Fig. |l} After reconstruction, the MDD 
has thrown out a ring of separate lumps of charge den- 
sity, with each lump containing one electron and having 
a radius somewhat larger than the magnetic length £b- 
This Chamon-Wen edge with broken rotational symme- 
try is located at a distance of ^ 2iB from the inner part 
of the droplet, which continues to form a (smaller) MDD 
with filling factor one in the dot center. With the recon- 
struction of the droplet, fractions of angular momenta 
between m — 8 and m — lA are moved to values above 
m = 19. The total angular momentum now is raised 
to M = 205.13. (Note that for the internally broken 
symmetry of the mean field solution, the total angular 
momentum can take fractional values.) 

The total current j(x, y) is plotted as a vector diagram 
in Fig. ||. It shows vortices which are located around the 
density maxima along the edge. 




FIG. 2. Real particle current j(r) for = 20 at 
B = 3.0T {left) and B = 3.5T {right), shown in the 
lower left quadrant of the charge density. The maxi- 
mum currents are about 0.028(nm ps)~^ (at 3.0T) and 
0.035(nm ps)~^ (at 3.5T). The shaded background in- 
dicates the charge density distribution: areas colored in 
magenta correspond to the regions with largest charge 
density 

For still higher fields, rotational symmetry is also bro- 
ken for the inner parts of the droplet (see lowest panel 
in Fig. |l|). The reconstructed density now forms a se- 
quence of rings, each consisting of well-separated max- 
ima in the electron density. Correspondingly, inside the 
droplet the real current ]{x,y) forms vortices (see right 
panel of Fig. ||) centered at the separate density max- 
ima along the rings. The localization of all electrons is 
associated with opening a large gap at the Fermi level. 
This had not been the case if only electrons along the 
edge were localized, as the inner part of the droplet then 
still was in the MDD phase. In the high-field limit, the 



"bumpy" electron density is consistent with the numeri- 
cal results of HF calculations by Miiller and Koonin pi[ . 
Performing exact calculations for < 6, Maksym |22] 
found localized states using a rotating frame. 




FIG. 3. As Fig. 0, but for A^ = 42 electrons and at a 
magnetic field B = 2.7T 

The broken symmetry of the Chamon-Wen edge is not 
limited to small dots. Figure || shows the density and 
angular momentum distribution of a dot with A^ = 42 
electrons. Again the edge consists of nearly localized 
electrons. The inner MDD at filling factor one is very sta- 
ble, with integer occupancy of angular momentum. By 
analyzing the Kohn-Sham single-particle orbitals we no- 
ticed that they fall into two discrete subsets, one forming 
the MDD and the other forming the broken-symmetry 
Chamon-Wen edge. This opens up a possibility for col- 
lective excitations localized at the edge. 




Magnetic Field B [T] 



FIG. 4. Different phases for a quantum dot with 
A^ = 6, 12, 16, 20, 30 and A^ = 42 electrons at a density 
roughly corresponding to = 2a|j. The colored areas 
indicate regions of magnetic fields where different phases 
(polarization, MDD, reconstruction and beyond) occur: 
The purple region indicates the MDD phase, while the 
magenta region corresponds to the Chamon-Wen edge 
regime. In the pink region localization extends over the 
whole droplet 

We finally study the formation of the MDD and its re- 
construction systematically as a function of the electron 
number. The resulting phase diagram is shown in Fig. ^. 
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We find three phase boundaries: the polarization tran- 
sition with the subsequent formation of the MDD, the 
broken-symmetry Chamon-Wen edge phase and finally 
the localization phase ||2l[ where the MDD completely 
disappears. The range of magnetic fields where the MDD 
is a stable ground-state becomes much more narrow with 
increasing N, as it has been noted in Refs. 0,^,^. 
The phase boundaries are in good agreement with re- 
cent experiments by Oosterkamp et al. j5|. (The fact 
that we find the transitions between different phases at 
smaller fields is due to a different electron density [Q). 
The experimental results also show systematic changes in 
the chemical potential after the first appearance of the 
Chamon-Wen edge. This might be related to the for- 
mation of rings and the stepwise disappearance of the 
MDD, which has a similar dependence of the number of 
electrons and the magnetic field. 

We have exclusively investigated the reconstruction 
in the regime where the droplet was completely polar- 
ized and refer to [ p5|j26| for a discussion about the ex- 
istence of spin textures along the edges. Performing 
fully unrestricted Hartree-Fock calculations, Karlhede 
and Lejnell |^5| as well as Franco and Brey found 
charge density wave-like modulations along the recon- 
structed edge of a Hall bar. These results are simi- 
lar to our finding of the broken-symmetry edge in finite 
droplets. 

In conclusion, we found that the maximum density 
droplet reconstructs into Chamon-Wen edge states with 
broken rotational symmetry in the internal coordinates. 
For larger fields reconstruction continues by sequential 
formation of rings up to an overall bumpy ground state 
density. The transitions with increasing B define dif- 
ferent phases of the electron liquid inside the droplet. 
We showed a phase diagram which displayed a system- 
atic separation between different phases of reconstruction 
both as a function of dot size and magnetic field. The 
trends for the shape of the phase boundaries between 
polarization transition, MDD, broken symmetry edge re- 
construction and finally localized states are rather 
clear and agree |2j] with the phase diagram recently mea- 
sured by Oosterkamp et al. [|]. 
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the TMR programme of the European Community un- 
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